Critical Current Behavior in Josephson Junctions with the Weak Ferromagnet PdNi 
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We have studied the variation of critical current in Superconductor/Ferromagnet/Superconductor 
(S/F/S) Josephson Junctions as a function of ferromagnet thickness using a weakly-ferromagnetic 
alloy Pd82Nii2. Measurements were performed for the thickness range 32 to 100 nm, over which the 
critical current density decreases by five orders of magnitude. The critical current density oscillates 
with a period of n x 4.1 ± 0.1 nm, and decays over a characteristic length of 8.0 ± 0.5 nm. There 
is no evidence of a crossover to a slower decay, which might indicate the presence of long-range 
spin-triplet pair correlations. We discuss possible reasons for their absence, including the possibility 
of strong spin flip scattering in PdNi. 

PACS numbers: 74.50.+r, 74.45.+C, 75.30.Gw, 74.20.Rp 
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I. INTRODUCTION 

When a superconducting (S) metal is placed in con- 
tact with a normal (N) metal, the properties of both 
metals are modified near the S/N interface. The result- 
ing superconducting proximity effect was widely stud- 
ied in the 1960's,- and again in the 1990's.-^ When the 
normal metal is replaced by a ferromagnetic (F) metal, 
the resulting physics is extremely rich, due to the very 
different order parameters in the two metals. There 
has been sustained interest in S/F systems over at least 
the past decade4>*i Many new phenomena have been ob- 
served, including but not limited to: oscillations in the 
critical temperature of S/F bilayers as a function of F 
layer thickness, 6 variations in the critical temperature 
of F/S/F trilayers as a function of the relative magne- 
tization direction of the two F layers^ and oscillations 
in the critical current of S/F/S Josephson junctions? 8 - 
The oscillatory behaviors are a direct result of the ex- 
change splitting of the spin-up and spin-down bands in 
the ferromagnet, which produce a momentum shift be- 
tween the up- and down-spin electrons of a Cooper pair 
that "leaks" from S into F. 9 In the clean limit, the spa- 
tial period of the oscillations is governed by the exchange 
length, £p = hvp/2E ex , where vf and E ex are the Fermi 
velocity and exchange energy of the F material, respec- 
tively. 

Unlike the S/N proximity effect, which extends to dis- 
tances of order 1 /im at sufficiently low temperature, 
the novel phenomena associated with the S/F proxim- 
ity effect persist only over very short distances - lim- 
ited cither by £p, by the mean free path, l e , or by 
their geometric mean, the dirty-limit exchange length 
%F = y/l^D / E ex ) , where D — vpl e /3 is the diffusion 
constant. These distance scales, characterizing the os- 
cillation and decay of superconducting correlations, tend 
to be extremely short in strong ferromagnets, and only 
moderately longer in the weakly-ferromagnetic alloys pre- 
ferred by some groups. For example, oscillations in the 
critical current in Nb/Co/Nb Josephson junctions have 
been observed with a period of 1.0 nm and a decay con- 
stant of 3.0 nm»i£ In Nb/Cu47Nis3/Nb alloy junctions, 



the observed period and decay constants were 11.0 nm 
and 1.3 nm, respectively^ 

In this context, it was very exciting when two theoret- 
ical groups predicted that a new form of superconduct- 
ing order, with spin-triplet pairing, could be induced by 
certain types of magnetic inhomogeneity in S/F systems 
consisting of conventional spin-singlet superconducting 
materials ! 11 ' 12 ^ 3 ' 14 Cooper pair correlations with spin- 
triplet symmetry are not subject to the exchange field 
of the ferromagnetic, since both electrons in the pair en- 
ter the same spin band in the ferromagnet. As a result, 
proximity effects due to such correlations should persist 
over long distances in a ferromagnetic, limited either by 
the temperature or by spin-flip and spin-orbit scatter- 
ing. Furthermore, the spin-triplet correlations predicted 
to occur in S/F systems are not of same type discov- 
ered recently in materials such as SrRuCU^ In SrRu04, 
the spin-triplet Cooper pairs satisfy the Spin-Statistics 
Theorem by having odd orbital angular momentum. In 
contrast, the spin-triplet correlations predicted to occur 
in S/F systems have even orbital angular momentum, 
and satisfy the Spin-Statistics Theorem by virtue of be- 
ing odd in time or frequency^ 

Experimental confirmation of the presence of spin- 
triplet correlations is not easy. In retrospect, several 
theorists have suggested their role in old experiments 
performed on mesoscopic S/F hybrid samples, where the 
data were interpreted in terms of a long-range supercon- 
ducting proximity effect . 17 ' 18 ! 19 More recently, there has 
been one report of a Josephson current in S/F/S struc- 
tures using CrC>2 as the F material, 20 where the dis- 
tance between the S layers was very long (several hun- 
dred nm) , so that the conventional spin-singlet supercur- 
rent should be exponentially suppressed. In a different 
work, 21 phase-coherent oscillations were observed in the 
normal resistance of a Ho wire connected to two super- 
conducting electrodes, again separated by a distance too 
large to support spin-singlet correlations. In both cases, 
the data were interpreted as being due to spin-triplet su- 
perconducting correlations in the ferromagnetic material. 
While these pioneering experiments are highly sugges- 
tive and tantalizing, the first suffers from large sample- 
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to-sample fluctuations in the magnitude of the observed 
supercurrent^ while both lack direct evidence for spin- 
triplet correlations. 



The goal of this work is to study S/F/S Josephson 
junctions where the thickness of the F layer is increased 
systematically, from a thin regime where the supercur- 
rent is likely to be dominated by spin-singlet correla- 
tions, to a thick regime where the smaller but longer- 
range spin-triplet supercurrent takes over. We chose Nb 
as the S material because its large critical temperature 
allows us to make measurements at 4.2 K, and hence 
measure a large number of samples. For the F material, 
we chose Pdi-zNi,,, alloy with a Ni concentration of 12 
atomic %, a material studied extensively by the group of 
M. Aprili 2 ^ and otherSfZ^iSLSLZS but only with PdNi 
thicknesses less than 15 nm. Our choice of this particular 
weakly-ferromagnetic alloy was based on two considera- 
tions: 1) Using a weakly-ferromagnetic material allows 
us to increase the thickness of the F layer without in- 
troducing an overwhelming amount of intrinsic magnetic 
flux inside the junctions. (This issue will be discussed fur- 
ther in Section III below.) 2) Some weakly-ferromagnetic 
alloys suffer from strong spin-orbit and/or spin- flip scat- 
tering, which are likely to destroy both spin-singlet and 
spin-triplet correlations. For example, Ryazanov and 
co-worker a 8 i 29 have found that the Josephson current in 
S/F/S junctions using Cui-^Ni^ alloy with x = 53 at.% 
decreases exponentially over a length scale of only 1.3 
nm. Moreover, that length scale is much shorter than 
the one characterizing the critical current oscillations, a 
fact that implicates strong spin-flip scattering in CuNi 
alloyi 9 ' 14 i 30 In contrast, while Kontos et al. 23 found that 
the critical current in S/F/S junctions made with PdNi 
alloy decays over a length scale only slightly longer, 2.8 
nm, they found that the oscillations and decay are gov- 
erned by the same length scale, which may imply that 
spin-flip and spin-orbit scattering are weak in this mate- 
rial. 



The paper is organized as follows. In Section II we 
discuss sample fabrication methods and characterization 
of the PdNi alloy. Section III discusses the characteriza- 
tion of our S/F/S Josephson junctions, with particular 
attention to the magnetic-field dependence of the criti- 
cal current (the so-called "Fraunhofer pattern"). Section 
IV presents the main results of the paper, namely the 
critical current vs. PdNi thickness of our S/F/S Joseph- 
son junctions. Section V discusses the various theoretical 
works on S/F/S junctions, and the physical parameters 
that one can extract from fitting theoretical formulas to 
the data. Our interpretation of the results is presented 
in Section VI. Finally, we conclude with suggestions for 
future directions. 



II. EXPERIMENTAL 
A. Sample Fabrication 

Substrates were silicon chips of dimension 12.7 x 
12.7 mm. Preparation for deposition was per- 
formed in a clcanroom, to minimize the presence 
of dust particles which could lead to shorts in the 
Josephson junctions. A multilayer consisting of 
Nb(150)/PdNi(d P(Wi )/ Nb (25)/Au(15) ( w j tn a n thick- 
nesses in nm) was deposited using magnetically-enhanced 
triode dc sputtering in an Ar plasma pressure of 2.5 
mTorr after obtaining a base pressure of 2 x 10" 8 Torr or 
better. A mechanical shadow mask was used to create the 
multilayer strip of size 0.16 x 10 mm 2 . The thin Au pro- 
tective layer prevents oxidation of the top Nb layer during 
further processing steps. The PdNi thickness, dpdNi was 
varied from 32.5 nm to 100 nm, typically in 5 nm steps. 
A subset of samples was fabricated with more closely- 
spaced thicknesses in the range from 58 nm to 75 nm, 
in order to demonstrate the minima in critical current 
indicative of the — n transitions. 

The multilayer was patterned using photolithography 
to create circular photoresist pillars with diameters of 
10, 20, 40, and 80 fim on each substrate. Care was 
taken to ensure the presence of undercut in the resist 
profile. Either a trilayer photolithography consisting of 
two photoresist layers separated by a thin metallic layer 
or a single layer photolithography using chlorobenzene 
yielded large reliable undercuts in the resist profile. The 
large variation in area allowed us to have a large dynamic 
range for the critical current measurements. In practice, 
however, the success rate of the 80 \xm pillars was very 
low, possibly due to the presence of dust particles dur- 
ing one of the fabrication steps taking place outside the 
cleanroom. 

The photoresist pillars acted as a mask to protect the 
multilayer below them while the rest was ion milled. The 
multilayer was milled down to the middle of the ferromag- 
netic layer, thus completely removing the top Nb layer 
yet not exposing the bottom Nb layer, to prevent the 
possibility of back-sputtered Nb depositing on the sides 
of the circular pillars. Nearly 200 nm of SiO x was then 
deposited to insulate the bottom Nb from the top Nb 
leads. Lift off of the photoresist pillars was then done us- 
ing Remover PG. This was followed by a slight ion milling 
of the top Au layer to ensure a clean interface. The top 
Nb lead of thickness 150 nm was then sputtered in the 
end. The thin Au layer becomes superconducting due to 
the proximity effect, as it is sandwiched between two Nb 
layers. A schematic of a complete Josephson junction is 
shown in Fig. 1. 

B. Characterization of PdNi Alloy 

The Ni concentration of our PdNi alloy was estimated 
by three different methods. Energy dispersive X-ray 
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FIG. 1: (color online). Schematic of a ferromagnetic Joseph- 
son junction. 



analysis (EDAX), performed on a 1.5/im thick PdNi film, 
yielded a Ni concentration of 12 ± 0.5%. (A thick film 
was used for this measurement to increase the signal- 
to-noise ratio for the Ni-K peak. Similar measurements 
performed on 200-nm thick PdNi film yielded a similar 
concentration value provided the signal was accumulated 
for long enough times.) To corroborate this value for the 
concentration, the magnetization M vs. temperature of 
a 100-nm thick PdNi film was measured using a Quan- 
tum Design SQUID magnetometer (see Fig. [2]). A clear 
change in the slope is seen around the Curie tempera- 
ture of about 175 K, independent of whether the mag- 
netic field is applied in-plane or out-of-plane. This Curie 
temperature corresponds to a Ni concentration of 12%, 
according to earlier work by BeiHe^ cited by Kontos<22 

Fig. © shows M vs. H at T = 10K for the same 100- 
nm thick PdNi film. The magnetization curve is more 
rounded, and has smaller remanent magnetization, when 
the field is applied in-plane, indicating that the mag- 
netic anisotropy of PdNi films is out-of-plane. Similar 
measurements on PdNi films of thickness 30 and 60 nm 
also indicate out-of-plane anisotropy, but somewhat less 
pronounced than in the 100-nm film. The out-of-plane 
anisotropy of PdNi surprised us initially, because the 
strong shape anisotropy of thin films usually dominates 
over magnetocrystalline anisotropy, resulting in overall 
in-plane anisotropy. Out-of-plane anisotropy has been 
observed in several materials, however, including CuNi 
alloyj 33 i 34 Recently, we learned that Aprili has also ob- 
served out-of-plane anisotropy in PdNi films. 35 The im- 
plications of PdNi's out-of-plane anisotropy on our work 
will be discussed further at the end of this paper. 

The saturation magnetization of PdNi measured out- 
of-plane at 10 K is M aat = 150 e mu /cm 3 = 0.23 /is/atom 



3lH32l |. this corresponds 



(see Fig. According to refs. 
to a Ni concentration of nearly 12%, in agreement with 
the determination from the EDAX measurements and 
Curie temperature. 



C. Ic measurement 

The normal-state resistances of our Josephson junc- 
tion pillars vary from 2.4 to 152 /zfi, depending 
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FIG. 2: (color online). Magnetization vs. temperature for a 
100-nm thick PdNi film, grown on 150 nm of Nb to have iden- 
tical crystalline properties as our Josephson junction samples. 
The film was first cooled in zero field to 10 K, then magne- 
tized by applying an in-plane field (black circles) or out-of- 
plane field (red squares) of 5 kG, then the magnetic moment 
was measured in field while heating to room temperature. 
The cooling curves (not shown) were also measured, and were 
found to follow the same curves as during heating except for 
occasional jumps seen for the out-of-plane field case. 



150 




-4000 -2000 2000 4000 

Field (Oe) 



FIG. 3: (color online) . Magnetization vs. in-plane field (black 
circles) and out-of-plane field (red squares) for the same 100- 
nm thick PdNi film shown in Fig. [2] measured at 10 K. 



mostly on the pillar area, and, to a lesser extent, on the 
PdNi thickness. Resistances in this range require an ex- 
tremely sensitive low noise measurement technique which 
is provided by using a superconducting quantum inter- 
ference device (SQUID) as a null detector in a current- 
comparator circuit.— All the four probe measurements 
were performed at 4.2 K by dipping the probe into a 
liquid Helium dewar equipped with a cryoperm shield. 
Each chip had Josephson junctions of diameters 10, 20, 
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40, and 80 /im. All measurements reported here were per- 
formed on junctions having Josephson penetration depth, 
-V/ = ($0/(2^^0 J c (c?f + 2Al)) 1 / 2 larger than one-quarter 
of the junction diameter w. (^o — h/2e is the supercon- 
ducting flux quantum, J c is the critical current density, 
and Al is the London penetration depth, equal to about 
86 nm in our sputtered Nb.) This ensures uniform cur- 
rent density in the Josephson junction* If Aj <C w, 
then the flux is screened from the center of the junction 
(Meissner effect) and the effect of the self field cannot be 
neglected. 

III. CHARACTERIZATION OF JOSEPHSON 
JUNCTIONS 

Fig. Q] shows a I-V relation typical for our S/F/S 
Josephson junctions. The curve follows the standard 
form for large-area, overdamped junctions: 

y(J) = ^Re[(/ 2 -J c 2 ) 1/2 ] (1) 

Occasionally, we find that the I-V curves are shifted hor- 
izontally, so that the critical current is not exactly the 
same in the positive and negative current directions. In 
such cases, we average the critical currents in the two 
current directions. 




Current (mA) 

FIG. 4: (color online). Voltage vs. current for a Josephson 
junction with diameter 20 /im and dpdNi= 62 nm and H ext 
— -19 Oe. The red solid line is a fit to Eqn. Q] 

One of the best ways to characterize Josephson junc- 
tions is to observe the modulation of critical current as 
a function of magnetic field H ex t applied perpendicular 
to the current flow direction in the junction. In nonmag- 
netic square junctions, the pattern so obtained is called 
the Fraunhofcr pattern, due to its similarity to the pat- 
tern produced in single-slit diffraction of light. Obser- 
vation of a good Fraunhofer pattern demonstrates that 
the supercurrent is uniform across the junction area, and 



that there are no short circuits in the surrounding SiO x 
insulator. 

In nonmagnetic Josephson junctions with circular cross 
section and negligible screening (Aj > w), the magnetic- 
field dependence of the critical current is given by 

2 x Ji(^) 
J($)=J C (0) \«° ; , (2) 

where 7 C (0) is the critical current in the absence of mag- 
netic field, Ji is the Bessel function of the first kind of 
order 1, and <f> = H ext (2XL + d)w is the magnetic flux 
penetrating the middle of the Josephson junction, with 
A l the London penetration depth, w the diameter of the 
circular junction and d the thickness of the barrier. This 
pattern is called an "Airy pattern." The pattern is quali- 
tatively similar to the Fraunhofcr diffraction pattern, but 
the first minima are spaced more widely apart in field 
than the subsequent minima. 

In samples containing a ferromagnetic barrier, one 
must include the intrinsic flux of the magnetic barrier 
on the Fraunhofer pattern. If the magnetization M is 
uniform throughout the junction, the intrinsic magnetic 
flux is given by $j? = fioMdpw, with dp the F-layer 
thickness and w the cross section width (equal to the 
diameter for circular junctions). In that case, the total 
magnetic flux through the F-layer is given by 

<S>tot = noMd F w + H ext (2\ L + d F )w (3) 

In macroscopic samples the magnetization breaks into 
domains, and Eqn. ([3]) is not valid. Instead, one must 
integrate the current density across the area of the junc- 
tion, taking into account the spatial dependence of the 
magnetic vector potential A due to the domains. In the 
Coulomb gauge, one takes A pointing along the current 
direction (z). The gauge- invariant phase difference across 
the junction includes a term proportional to the line in- 
tegral of A from deep inside one superconducting contact 
to deep inside the other^ The resulting expression for 
the supercurrent is: 

la = J J(x, y)dxdy (4) 

where 

J{x,y) = J c sin(a + ^- J A(x,y) .dl) (5) 

The term containing the vector potential performs a ran- 
dom walk as one moves across the sample, due to the 
domains pointing in random directions. If the magnetic 
domains are very small and/or the magnetization is very 
weak, then the vector potential term stays near zero 
in all parts of the junction, and the critical current is 
hardly affected. If, however, the magnetic domains are 
large and/or have large magnetization as in the case of 
strong ferromagnet, the contribution to the phase due to 



5 



the vector potential deviates far from zero as it crosses 
even a single domain, thus severely suppressing the crit- 
ical current. This can lead to complete destruction of 
the Fraunhofer pattern. This is clearly seen in Fig. [5j 
which shows data for an S/F/S junction of diameter 10 
fim, with an 11-nm thick Ni layer. (Similar random- 
looking "Fraunhofer patterns" have also been seen by 
other groups studying S/F/S Josephson junctions^) In 
principle, a regular Fraunhofer pattern can be recovered 
if the sample is completely magnetized, by applying a 
magnetic field in the opposite direction to the magneti- 
zation such that the total flux given by Eqn. §3§ is zero. 
In that case, one should observe a regular Fraunhofer 
pattern shifted in field by an amount: 




FIG. 5: (color online). Critical current vs. in- plane magnetic 
field for a Nb/Ni/Nb circular Josephson junction of diameter 
10 /im, with djvi = llnm. The black points (squares) were 
measured in the virgin state, whereas the red points (circles) 
were measured after magnetizing the sample in an external 
field of +1 kOe. The random pattern arises due to the intrin- 
sic magnetic flux of the complex domain structure of the Ni 
layer. 

The above argument holds only if the coercive field of 
the magnetic material is large enough so that the magne- 
tization stays nearly uniform even in the presence of the 
applied field, H s hift, pointing in the opposite direction. 
For the case of the Ni sample shown in Fig. the largest 
peak in the critical current vs. field after magnetization 
is found near 160 Oe, whereas the expected shift calcu- 
lated from the known saturation magnetization of Ni is 
about 207 Oe. The discrepancy is caused by some rota- 
tion of the magnetization in the domains or some domain 
wall motion as H ext approaches the coercive field, which 
we measured to be approximately 180 Oe in a 9-nm Ni 
film. 



To avoid the distortion of the Fraunhofer pattern, there 
are several options for the study of ferromagnetic Joseph- 
son junctions: 1) Use ultrathin samples: Under this con- 
dition the flux enclosed in the junction due to a single 
magnetic domain is much less than one flux quantum, 
and one can safely ignore the contribution to the flux 
from the magnetization. This option is not available to 
us, because our goal requires us to work with thick fer- 
romagnetic layers. In addition, thin magnetic films of- 
ten have magnetically "dead" layers on each side, which 
pose additional problems for ultra thin samples; 2) Use 
samples with ultra-small lateral dimensions to reduce the 
contribution to the total magnetic flux from the magneti- 
zation, and to control the domain structure: This method 
has been pursued by Blamire and co-workers^ and also 
by Strunk and co-workersj« using strong ferromagnets. 
The disadvantage of this approach is that it becomes less 
effective as the thickness of the ferromagnetic layer is 
increased. In addition, it restricts the possibilities to 
introduce magnetic inhomogeneities that are naturally 
present due to the domain structure in devices of larger 
dimension, and which may be crucial for inducing the 
predicted spin-triplet superconducting correlations dis- 
cussed in the Introduction; 3) Work with materials that 
either have weak magnetization, small domain size, or 
both: As discussed previously, this approach has been 
used by Ryazanov and co-worker o 8 ' 29 who worked with 
CuNi alloy, and by Aprili and co-workers ^ who worked 
with PdNi alloy. In CuNi alloy, Ryazanov demonstrated 
that the magnetization of the CuNi makes very little con- 
tribution to the total magnetic flux in his samples^ 2 - He 
did this by comparing the Fraunhofer pattern for a de- 
magnetized sample with the pattern for the same sample 
uniformly magnetized. The latter pattern was shifted by 
a constant field, while the maximum value of the critical 
current was nearly unchanged. That shows that the in- 
tegrated vector potential stayed close to zero everywhere 
in the sample; 4) Engineer the F layer to have zero net 
magnetic flux, for example, by using a "synthetic antifer- 
romagnet." We are currently exploring this option, and 
will report it in a future publication^ 

For this work, we have chosen the third option, but 
with PdNi alloy rather than CuNi alloy as our weak fer- 
romagnet. As discussed in the introduction, there is evi- 
dence of strong spin-flip scattering in CuNi alloy, whereas 
the situation in PdNi alloy is less clear. It should be em- 
phasized that the magnetism in PdNi is quite different 
from the magnetism in CuNi. Because Pd is nearly ferro- 
magnetic itself according to the Stoner criterion, it takes 
only a small concentration of Ni to make the alloy ferro- 
magnetic. One might expect then that the magnetism is 
more uniform in PdNi than in CuNi, where magnetism 
results from large Ni clusters. 

As a final note regarding Nb/PdNi/Nb junctions, we 
found that using PdNi layers thinner than 30nm resulted 
in very large critical current densities - so large that 
even our smallest pillars (w — lOfim) were in the regime 
A,/ <C w. Kontos et al£^ circumvented that problem by 
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introducing an additional insulating layer in their junc- 
tions to reduce J c and hence increase Aj. Because our 
interest is in studying junctions with large dp, we omit- 
ted the insulating layer. This choice limited our study to 
junctions with dp > 30nm, which have small enough J c 
so that A j > w/4 for our smallest diameter pillars. 
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FIG. 6: Critical current vs. applied magnetic field (Fraun- 
hofer patterns) obtained for Nb/PdNi/Nb Josephson junc- 
tions with different thickness of PdNi interlayer a) 35nm, b) 
50nm, c) 70nm, d) 85nm. The pillar diameters w are 10, 10, 
f 0, and 20 fim, respectively. 



IV. S/F/S JOSEPHSON JUNCTIONS WITH 
PDNI: RESULTS 
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FIG. 7: (color online). Critical current density vs. dpdNi 
for all of our Nb/PdNi/Nb Josephson junctions. Each data 
point represents the average over multiple pillars on the same 
substrate, with the error bars as the standard deviation of the 
mean. Black points (squares) are virgin state data, while red 
points (circles) were measured after magnetizing the samples 
(see Fig. [TO}. The solid line is a fit of Eqn. © to the 
virgin state data, while ignoring the last two data points with 
dpdNi = 95 and 100 nm. 



neighborhood of 35, 48, 60, and 75 nm. Figs. [5]and[l 
show J c vs G?i? on linear axes, where the local minima 
are more clear. Such local minima have been observed 
in S/F/S junctions made with several different ferromag- 
netic materials^ i 10 ' 23 i 44 ' 45 i 46 ' 47 i 48 and signify transitions 



between standard junctions and 7r-junctions. 



Fig. [S]shows I c vs. H ext data for Nb/PdNi/Nb Joseph- 
son junctions with 35, 50, 70, and 85 nm of PdNi. The 
critical current has minima whenever an integer number 
of flux quanta penetrate the junction. The extremely low 
values of critical current at the minima indicate the ab- 
sence of any shorts in the insulating material surrounding 
the Josephson junctions. Similar measurements were per- 
formed on samples for which the PdNi thickness varied 
from 35 nm to 100 nm. The maximum current density, 
J c , obtained for all such devices is plotted vs. thickness 
in Fig. [7] This figure represents the main result of this 
work. The error bars represent the standard deviation 
of the mean of the results obtained from several pillars 
on the same substrate. (For the smaller values of dp, 
we measured primarily the pillars of diameter 10 and 20 
/im, whereas for the larger values of dp, we measured the 
20 and 40 pillars.) The critical current density de- 
creases exponentially over five orders of magnitude over 
this range of PdNi thickness. To our knowledge, these 
data represent the widest range of ferromagnet thickness 
in S/F/S Josephson junctions studied to date. 

Fig. [7] shows that J c does not decrease monotonically 
with dp, but rather exhibits local minima with dp in the 
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FIG. 8: (color online). Linear plot of J c vs. dpdNi for the 
thickness range 32.5 to 58 nm. The line is the same fit shown 
in Fig. [3 

Measurements of I c vs. H ext were also performed on 
the samples after magnetizing them by applying an in- 
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FIG. 9: (color online). Linear plot of J c vs. dpdNi for the 
thickness range 55 to 85 nm. The line is the same fit shown 
in Fig. 

plane field of 5 kOe. The resulting Fraunhofer patterns 
(see Fig. ITU]) are shifted in field to a point where the flux 
due to the external field cancels out the flux due to the 
intrinsic magnetization. The maximum critical currents 
obtained for magnetized samples match closely with the 
virgin state data. This indicates that the domains in 
PdNi alloy are relatively small, so that the total excursion 
of the integrated vector potential as one crosses a domain 
is much less than one flux quantum. 
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FIG. 10: Fraunhofer pattern in the virgin state (right) and 
after magnetizing (left) a Nb/PdNi/Nb Josephson junction 
with diameter w = 10p,m and dpdNi=47-5 nm. 

Fig. QTJ shows the field shift of the Fraunhofer pat- 
tern of the magnetized samples, vs. PdNi thickness cIf- 
As d,F increases, the field shift saturates at a value near 
200 Oe. As the PdNi thickness increases, there is an 



increasing tendency for the magnetization to rotate out 
of the plane, thereby decreasing its in-plane component. 
The solid line is a fit to the data of Eqn. [51 with the 
only free parameter being the remanent magnetization 
M = 54 emu/cm 3 . This compares with values of 76 and 
62 emu/cm 3 measured directly on PdNi films of thick- 
ness 30 and 60 nm, respectively. The red stars in Fig. 1111 
show the values of H s hijt calculated from Eq. using 
the values of M rem measured directly on PdNi films of 
thickness 30, 60, and 100 nm. The agreement with the 
field shifts of the Fraunhofer patterns is reasonable. 
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FIG. 11: (color online). Magnetic field shift of Fraunhofer 
patterns of magnetized samples vs. PdNi thickness d,F- The 
line is a fit of Eqn. ([6]) to the data points with (If < 72 
nm, with \l = 86 nm. The fit provides an estimate of 54 
emu/cm 3 for the remnant magnetization M rem of PdNi alloy. 
The red stars indicate the calculated H s hift using values of 
M rem measured in a SQUID magnetometer on PdNi films of 
thickness 30, 60, and 100 nm. (Fig. [3] shows the M vs. H 
data for the 100-nm film.) 

We have also measured the normal-state resistance of 
our samples at currents much larger than I c . A plot of the 
specific resistance ARn (area times resistance) vs. dpdNi 
is shown in Fig. 1121 The interface and bulk contributions 
to ARjy are given by: 

AR N = 2AR B + ppdmdp (7) 

where ppdNi is the resistivity of PdNi, dp is the thickness 
of the PdNi layer, and Rb is the Nb/PdNi boundary re- 
sistance. A linear fit to all of the data gives a boundary 
resistance of ARp = 1.95 ± 0.14/S!m 2 and a resistivity 
of PdNi, ppdNi = 82 ± 4nilm. There is some indication 
in the data that the slope increases for dp greater than 
about 75 nm. If we fit only the data for dp < 75 nm, then 
we find AR B = 2.39 ± 0.14f2m 2 and p PdNi = 62±5nfim. 
Independent measurements of the in-plane PdNi resis- 
tivity were performed on 200-nm thick films, using the 
van der Pauw method. Those measurements produced 
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the value ppdNi = 116 ± 2n£lm. It is plausible that the 
in-plane resistivity is larger than the perpendicular re- 
sistivity if the PdNi films grow in a columnar fashion, 
although measurements on other sputtered metals often 
find quite close agreement between these two measure- 
ment methods.— 

13.5- 1 
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FIG. 12: Area times normal-state resistance vs. dpdNi for 
all of our Josephson junction samples. The slope provides 
the resistivity of PdNi and the y-intercept provides twice the 
Nb/PdNi boundary resistance 



V. THEORY OF S/F/S JOSEPHSON 
JUNCTIONS 

There have been a large number of theoretical works 
dealing with S/F/S Josephson junctions. There are three 
energy scales whose relative size determines three distinct 
regimes. The energy scales are the exchange energy in the 
ferromagnet, E ex , the gap in the superconductor, A, and 
h/r, the inverse of the mean free time between collisions 
of an electron propagating in the ferromagnet. In all 
of the experimental work on S/F/S Josephson junctions 
published to date, including this work, E ex 3> A. There 
is a wide variation, however, in the size of H/t relative to 
those two energies. The true clean limit is expressed by 
E ex T 3> h, which also implies At 3> h. The intermediate 
limit is where E ex r 3> h but At -C h, whereas the true 
dirty limit is where both At <C H and E ex r <C H. These 
three regimes can also be characterized by the relative 
sizes of the three length scales: the electron mean free 
path, l e = vpT, the superconducting coherence length, 
£s = HDs/A, and the clean-limit exchange length dis- 
cussed earlier, £p = Hvf /2E ex . 

The simplest limit is the dirty limit, with the addi- 
tional constraint that the ferromagnetism is weak enough 
so that one can treat the spin-up and spin-down bands 
identically, i.e. as having the same Fermi velocity and 



mean free path. In this limit, the Usadel equation is 
valid. An expression for the critical current I c as a func- 
tion of ferromagnetic layer thickness cLf was first derived 
by Buzdin et al£Q In this regime, the oscillation and de- 
cay of I c as a function of cLf are both governed by a 
single length scale - the " dirty- limit" exchange length, 
£p = (hDp/Eex) 1 ^ 2 . Once dp exceeds the thickness 
dependence of I c takes the simple form: 

I c {d F ) = I c0 * exp(^) * | sin(^ + tt/4)| (8) 

?F ?/ 

In the presence of spin-flip or spin-orbit scattering, Eqn 
© is modified, and the length scales governing the decay 
and the oscillation are no longer equal^^^i The more 
general form can be written as 

I c {d F ) = I c0 *exp(^)* | sin(-^ + <f>)\ (9) 

?F1 ?F2 

In general, the effect of spin-flip or spin-orbit scattering 
is to shorten the decay length scale, £pi, relative to 
and to lengthen the oscillation length scale, £f2- In the 
presence of sufficiently strong spin-orbit (but not spin- 
flip) scattering, the oscillations disappear completely. An 
equation similar to Eqn. ((9]) has successfully been used 
to fit I c vs. dp data from S/F/S junctions containing 
CuNi alloy, with £pi = 1.3nm, and £f2 = 3.5nm?£ The 
very short value of £fi compared to £f2 was interpreted 
as implying that spin-flip scattering is strong in that ma- 
terial. 

Eqn. ([9]) can also be fit to our data, as shown in 
Figs. [7JH1 But in our case, the length scale governing 
the exponential decay (£fi = 8.0 ± 0.5nm) is consider- 
ably longer than the length scale governing the oscillation 
(£f2 — 4.1 ± O.lnm), hence the dirty- limit hypotheses 
that led to Eqn. ((9]) are not fulfilled. However, the con- 
dition £pi > £f2 has also observed in S/F/S Josephson 
junctions containing the strong fcrromagnets: Ni, Fe, Co, 
and NigoFe2o (also known as Permalloy, or Py).— Those 
materials have very large exchange energy, hence they 
are in the intermediate limit, with E ex r ^> H but still 
At <C h. Our observation of £fi > £f2 in S/F/S junc- 
tions with PdNi alloy suggest that PdNi may also be in 
the intermediate limit. 

The intermediate limit has been studied in several the- 
oretical papers. Bergeret et al^ provide a general for- 
mula for the critical current, valid both in the dirty limit 
and intermediate limit. In the intermediate limit the for- 
mula simplifies when the F-layer thickness is larger than 
the mean free path. Eqn. (20) from ref. [ID is: 

I c (d F ) (x irT \ exp — — 1+2wt ) 

Z^ A +U} dp / £p l e 

(10) 

where the sum is over the positive Matsubara frequencies, 
w m = 7rfcsT(2m + 1) with T the temperature. Asymp- 
totically at large values of dp, Eqn. (fTUjl is quite similar 
to Eqn. ([5]), with £pi = l e and £i?2 = £f- Nevertheless, 



9 



we have fit Eqn. (fTUj) directly to our data, and obtained 
the parameters £p = 4.0 ± O.lnm and l e = 11.2 ± lnm. 
Not surprisingly, the former agrees closely with the value 
of £f2 obtained from the fit of Eqn. © , while the value 
of l e obtained from Eqn. (|10[) is somewhat larger than 
the value of £fi obtained from the fit of Eqn. ((9J - the 
difference undoubtedly due to the sum over Matsubara 
frequencies in Eqn. (| 10[) . We will use the larger value of 
l e in the following discussion. 

An alternative theoretical approach was taken by 
Kashuba, Blanter, and Fal'ko, 52 who considered a model 
of S/F/S junctions taking explicit account of spin- 
dependent and spin- flip scattering in the F layer. In 
the intermediate limit, these authors find a result that 
is nearly identical to Eqn. IjlOp. 



VI. DISCUSSION OF RESULTS 

A. Estimate of mean free path and exchange 
energy in PdNi 

In many metals, it is straightforward to estimate the 
mean free path directly from the measured resistivity, 
as the product pl e is inversely proportional to the Fermi 
surface area, and has been tabulated for a large number 
of metals^ In PdNi alloy, however, the pl e product is 
not known. Some workers have tried to estimate pl e in 
PdNi from its value in Pd, but even that estimation is 
not straightforward, due to the complex band structure 
of Pd. To illustrate the difficulty, previous workers have 
quoted values of the pl e product as small as 
(refs. HHi) and as large as 4.0/fJm 2 (ref. M). The for- 
mer value is certainly too small, because it was calculated 
using the Einstein relation, a — n(Ep)e 2 D, with the den- 
sity of states at the Fermi level, n(Ep), obtained from the 
electronic specific heat coefficient 7 = (n 2 /3)k 2 3 n(Ep). 
The problem is that, in Pd, the specific heat is dom- 
inated by heavy holes on the open "jungle gym" por- 
tion of the Fermi surface, while the electronic transport 
is dominated by much lighter electrons on a part of the 
surface centered at the T-point. (It is common to refer to 
these electrons as "s-like"; but that is incorrect because 
they are strongly hybridized with the d bands. 54 ) Discus- 
sion of the Fermi surface in Pd has been given in several 
papersi 5 ^ 5 -^! 5 ^ Pinski et al£l state that the T-centered 
sheet of the Fermi surface carries the vast majority of the 
transport current - up to 97% at 10 K. 

We estimate the pl e product for PdNi in two ways. 
In his Ph.D. thesis* 3 - 2 - Kontos found that the resistivity 
of thin PdNi films with x « 12% varied linearly with 
inverse thickness 1/d once the films were thinner than 
about 8 nm, indicating that the mean free path is lim- 
ited by the film thickness. The slope of the graph gives 
the product pd = 1.7 f Vim 2 . We expect the pl e product 
to be within a factor of two of this value. The second 
method relies on the statement by Pinksi et al£^ that 
electronic transport in Pd is dominated by the electrons 



on the T-centered sheet. The total number of carriers 
in that band, the density of states at the Fermi level, 
the effective mass, and the Fermi velocity on that sheet 
have all been tabulated by Dye et al^ based on deHaas 
van-Alphen measurements of the Fermi surface of Pd. 
The values of those four quantities are: n = 0.375 car- 
riers/atom = 2.54 ■ 10 28 7tj- 3 ; n(E F ) = 0.189 states/(eV 
atom spin) = 2.56T0 28 eV r ~ 1 m -3 ; m* = 2.0m e ; and vf = 
0.6 ■ (27r/a) ■ {h/2m) = 5.6 • 10 5 m/s^ Using the Drude 
formula, a — ne 2 T/m* , one finds pl e = 1.55/f2m 2 , while 
using the Einstein relation one finds pl e = 1.31/f2m 2 . 
(The slight difference between the two values is likely due 
to the non-parabolic character of the T-centered sheet.) 
These values are close to the value pl e — 1.7 f Vim 2 esti- 
mated from the thickness dependence found by Kontos. 
If we take our own measured resistivity, p — 82nf2m, and 
the value of the mean free path from the J c vs. dp fit, 
l e = 11.2nm, we obtain pl e = l.OfSlm 2 , which is not too 
far from the estimates given above. (The only mystery 
we can not explain from this analysis is the very low value 
of the Fermi velocity, vp = 2.0-10 5 m/s, measured by Du- 
moulin et al^ in a proximity effect experiment between 
Pb and Pd. It is unclear why that experiment measures 
the low Fermi velocity of the open hole sheet rather than 
the higher Fermi velocity of the T-centered sheet.) 



From the period of oscillation of J c vs. dp, we found 
£.F2 = 4.0 ± O.lnm. Using the Fermi velocity of the 
dominant carriers, vp — 5.6 • 10 5 m/s, gives an estimate 
for the exchange energy in our PdssNi]^ alloy of E ex — 
Hvp/2^p2 = 44meV. This value is somewhat higher than 
values quoted by previous worker a 23 i 24 i 25 i 26 i 27 , but those 
earlier estimates were either made using the smaller value 
ofvF, or using the diffusive formula £ F — \J~(KD j E ex ). A 
more meaningful comparison is of the length scales t^p\ 
and £f2 found in different experiments. For example, 
Kontos et al^ found £fi rj £,f2 = 2.8nm in their study 
of S/I/F/S Josephson junctions with a PdNi alloy of sim- 
ilar concentration to ours. The values of £f2 in their 
experiment and ours are rather close to each other, but 
the values of £fi are not. It is not clear if that discrep- 
ancy is significant or not. The thickness range covered 
in the earlier work was 4.5 - 14 nm, whereas the range 
we covered was 32.5 - 100 nm. If the mean free path in 
the PdNi alloy is indeed in the range of 11 nm, then the 
samples studied by Kontos et al. were in the crossover 
regime with dp m l e , where the thickness dependence 
has not yet obtained the asymptotic exponential decay, 
J c oc exp(— dp/l e ). But then one would expect a less 
steep decay of J c with dF, rather than a more steep de- 
cay. Perhaps a more relevant observation is simply that 
the PdNi films deposited in different laboratories may 
have different polycrystalline structures, and hence very 
different mean free paths. A summary of the parameters 
estimated by previous workers, as well as by our work, is 
given in Table I. 
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Source 


Experiment 


Ni cone. 


VF 


pie 


P 




fa 


£F2 


Eex 


formula to 


Tcurie 






(at. %) 


(10 5 m/s) 


(fflm 2 ) 


(nfim) 


(run) 


(nm) 


(nm) 


(meV) 


extract E ex 


(K) 


22 


S/F DOS 


10 


2.0 


4 


- 


~ cLf 


2.8 


2.8 


35 


(hD/E ex )^ 2 


100 


23 


S/I/F/S 


12 


2.0 


4 


- 


~ cLf 


3.5 


3.5 


13 


{hD/E^) 1 ' 2 


260 


25 


T c of S/F 


14 


2.0 


4 


240 


16.6 


6 


6 


15 


(hD/E ex )^ 2 


156 


26,27 


T c of S/F 


14 


2.0 


0.96 


240 


4 


3.4 


3.4 


13 


(hD/E ex )^ 2 


185 


24 


T c of S/F 


15 


2.0 


0.3256 


220 


1.48 


3.5 


3.5 


13 


{hD/E^) 1 ' 2 


260 


28 


S/F/S 


18 


2.0 








2.4 


2.4 


52 


(hD/E ex )^ 2 


200 


this work 


S/F/S 


12 


5.6 


1.0 


82 


11 


8.0 


4.0 


44 


hv F 1 E ex 


175 



TABLE I: PdNi parameters from several groups. The superconductor used in all the above experiments was Niobium. Note 
that most workers have used the diffusive formula, £f2 = {hD / E ex ) x ^ 2 , to extract E ex from the measured value of £f2, whereas 
we have used the ballistic formula, £f2 = Hvf /2E ex . Our choice of Fermi velocity, vf, is discussed in the text and in ref. l58l . 



B. Spin-triplet superconducting correlations? 

One of the primary goals of this work was to search for 
signs of spin-triplet superconducting correlations in our 
samples. At first glance, the data in Fig. ([7]) show no 
sign of spin-triplet superconducting correlations, which 
might manifest themselves as a crossover to a less-steep 
exponential decay of J c at large values of dp. It is in- 
triguing, however, that the length scale characterizing the 
exponential decay of J c in our samples, £fi = 8 — llnm, 
is substantially longer than that observed previously in 
shorter S/F/S junctions with PdNi^ Could it be that 
we are already observing the triplet Josephson effect 
throughout the whole range of dp reported here? The 
strongest evidence against such an interpretation is the 
nearly-periodic set of local minima we observe in J c vs. 
dp, shown in Figures ((7][9|). We believe that those lo- 
cal minima signal crossovers between 0-junctions and ir- 
junctions, which are due to the effect of the exchange 
splitting on spin-singlet superconducting correlations. A 
Josephson supercurrent dominated by spin-triplet corre- 
lations would not exhibit such minima, but rather would 
decay monotonically with increasing dp. Nevertheless, to 
rule out the spin-triplet hypothesis definitively would re- 
quire stronger evidence that the local minima we observe 
truly represent — n crossovers, rather than an unlucky 
distribution of uncertainties in the data that mimics a 
periodic set of local minima. 

There are three ways one could make a more stringent 
test that the local minima observed in our J c vs. dp data 
are due to — it crossovers: 1) A direct measurement of 
the current-phase relationship of the junction^ 2) Ex- 
tension of J c vs. dp measurements to smaller values of 
dp, to see if the slope of the J c vs. dp semi- log plot 
changes to a value close to that measured by Kontos et 
al. This would require reducing the lateral size w of our 
junctions, so as to maintain the condition Aj > u>/4; 3) 
Measurement of J c vs. temperature T for samples very 
close to a — it crossover. In S/F/S junctions with very 
weak ferromagnets, the — n crossover has been observed 
in the T-dependence of J e 4 i 29 i 45 i 48 As the exchange en- 
ergy increases, however, the thickness range over which 



one can see a non-monotonic T-dependence of J c gets 
progressively narrower 4£ Each of these checks presents 
its own set of challenges, and represents a possible direc- 
tion for future work. 



C. Spin- flip scattering in PdNi 

Assuming that our data are not the result of spin- 
triplet correlations, we would like to know why not. 
There are several factors that may contribute. First and 
foremost, strong spin- flip scattering, if it exists in our 
PdNi alloy, would destroy the triplet correlations. To 
address this issue, we have independently tried to mea- 
sure directly the spin diffusion length, l s f, in PdNi alloy 
using techniques borrowed from the Giant Magnetore- 
sistance (GMR) community. A complete discussion of 
those measurements is given elsewhere^ Here we note 
the most salient results. First, the value of the spin diffu- 
sion length obtained, l s f = 2.8nm, is extremely surpris- 
ing given the much longer length scale characterizing the 
decay of the Josephson supercurrent in the present work. 
Normally, one assumes that spin-flip and spin-orbit scat- 
tering processes occur on length scales much longer than 
the mean free path, justifying the diffusive model used 
in discussing the spin diffusion length^ Hence a mea- 
surement of a spin diffusion length several times shorter 
than the mean free path is difficult to interpret. We be- 
lieve that the very short value of l s f observed using a 
GMR spin valve^ may be due to the out-of-plane mag- 
netic anisotropy of PdNi discussed earlier. The analysis 
of the spin-valve data of ref. [6l| assume that the mag- 
netizations of the two magnetic layers are either parallel 
or antiparallel to each other. Our recent discovery of the 
out-of-plane magnetic anisotropy in PdNi casts doubt on 
this assumption. It is possible that the magnetization of 
PdNi is inhomogeneous on very short length scales, which 
may then destroy spin memory by rotating the spin on a 
length scale smaller than the mean free path. Clarifying 
this issue will require further experiments. 

A second possible reason we do not observe signs 
of spin-triplet superconducting correlations is that the 
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length scale characterizing the magnetic inhomogeneity 
in PdNi might not be comparable to the Cooper pair 
coherence length £ s in Nb. Let us refer to the length 
scale characterizing the magnetic inhomogeneity as S m . 
If S m <C £ s , then a Cooper pair will experience the mag- 
netization averaged over the length £ s , and triplet corre- 
lations will not be produced efficiently. If 5 m ;§> £ s , then 
a typical Cooper pair experiences no magnetic inhomo- 
geneity. The coherence length in our Nb is known to be 
about 13 nm.— Unfortunately, the typical domain sizes 
and domain wall widths in our PdNi alloy are not known. 
Because the Curie temperature of PdNi is well below 
room temperature, obtaining information about mag- 
netic structure requires a low temperature magnetic vi- 
sualization technique, such as low-temperature magnetic 
force microscopy (MFM) or Bitter decoration. And even 
then, the former method is not well-suited to weak ferro- 
magnets, because the magnetization of the MFM tip may 
influence the domain structure of the sample. Very re- 
cently the domain structure of CuNi alloy was measured 
using the Bitter decoration technique.— Those measure- 
ments confirmed the out-of-plane magnetic anisotropy of 
that material, and found a typical domain size of lOOnm 
in the virgin state or at the coercive field. Unfortunately, 
no such measurements have been performed on PdNi al- 
loy, to our knowledge. Clearly, a thorough study of the 
magnetic domain structure of PdNi would help clarify 
this issue. 



VII. FUTURE DIRECTIONS 

The most urgent work needed in the future is a strong 
verification (or repudiation) that the local minima in our 
data do indeed represent — ir crossovers, rather than 
sample-to-sample fluctuations. Looking further ahead, 
to have any hope of nailing down the elusive spin-triplet 
supercurrent in S/F/S junctions will require better char- 
acterization of magnetic materials. On the one hand, the 
spin diffusion length is a crucial parameter, as it lim- 



its the spatial extent of spin-triplet correlations. For- 
tunately, the spin-diffusion length has been measured 
in some ferromagnetic materials,— but more work is 
needed. Of equal important is information about the 
typical length scales characterizing the domain structure 
of ferromagnetic thin films. This is a complex issue, as ei- 
ther the domain size or domain wall width may be impor- 
tant. For example, in a film where neighboring domains 
have antiparallel magnetization, the long-range triplet 
component is generated only in the domain walls~ 4 hence 
it is crucial that the domain wall width be comparable to 
£ s . If, however, neighboring domains have non-collincar 
magnetizations, then it would seem that the long-range 
triplet could be produced even if the domain walls are 
very thin, as long as the typical domain size is compara- 
ble to £ s . Unfortunately, domain sizes in polycrystalline 
films are not known a priori^ 

VIII. CONCLUSIONS 

We have measured the critical current of Nb/PdNi/Nb 
Josephson junctions with PdNi thicknesses ranging from 
32.5 to 100 nm. The critical current drops by five orders 
of magnitude over this thickness range. The data ap- 
pear to have a periodic array of local minima, suggesting 
that the supercurrent is of the conventional spin-singlet 
type over the entire thickness range. We have discussed 
possible reasons for the absence of induced spin-triplet 
correlations, such as spin-flip scattering or a poor match 
of length scales between the magnetic domains and the 
superconducting coherence length. 
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